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Abstract. Recurrent novae occurring in symbiotic binaries are candidate 
sources of high energy photons, reaching GeV energies. Such emission is a 
consequence of particle acceleration leading to pion production. The shock 
between matter ejected by the white dwarf, undergoing a nova explosion, and 
the wind from the red giant companion is responsible for such a process, which 
mimics a supernova remnant but with much smaller energetic output and much 
shorter time scales. Inverse Compton can also be responsible for high energy 
emission. Recent examples are V407 Cyg, detected by Fermi, and RS Oph, 
which unfortunately exploded in 2006, before Fermi was launched. 

Key words: white dwarfs, novae, symbiotic stars, supernovae, X- and gamma- 
ray astronomy 

1. INTRODUCTION 

White dwarfs are the endpoints of stellar evolution of stars with masses smaller 
than about 10 Mq. In white dwarf stars there is not nuclear energy available 
anymore. Their chemical composition is either a CO or ONe mixture (it can also 
be pure He). Their typical masses are ~ 0.6 Mq and their maximum mass is the 
Chandrasekhar mass (~ 1.4 Mq). White dwarfs cool down to very low luminosities 
(typically L ~ 10~*'^ Lq) when isolated, whereas they can explode when they are 
accreting matter in close binary systems. 

There are two main types of binary systems where white dwarfs can accrete 
matter and subsequently explode. The most common case is the so-called cata- 
clysmic variable, where the companion is a main sequence star. In this system, the 
mass transfer occurs via Roche lobe overflow, typical orbital periods are ranging 
from hours to days, and orbital separations are of several 10^" cm. Hydrogen burn- 
ing in degenerate conditions on top of the white dwarf leads to a thermonuclear 
runaway and a (classical) nova explosion. Masses ^ 10~^-10~'' Mq are ejected 
at velocities of several lOO's or even lOOO's of km s~^. A nova explosion does 
not disrupt the white dwarf (as occurs in thermonuclear - or type la - supernova 
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Figure 1: Time-dependence of forward shock velocity in the 2006 outburst of RS 
Oph as deduced from FWZI of IR emission lines (filled squares and circles: Das et 
al. 2006, filled triangles: Evans et al. 2007) and X-ray measurements of the post 
shock temperature (open triangles: Sokoloski et al. 2006, open squares: Bode et al. 
2006). The IR data can be modeled by (dashed line) Vs{t) = 4300(i/^i)"" km s-\ 
where ti = 6 days and a„ = 0(— 0.5) for t < ti{t > ti). 



explosions); therefore, after enough mass is accreted again from the companion 
star, a new explosion will occur. The typical recurrence time is 10^-10^ years. In 
our Galaxy, about 35 classical nova explosions occur every year. 

Another scenario is the so-called symbiotic binary, where the white dwarf ac- 
cretes matter from the stellar wind of a red giant companion. Typical orbital 
periods for these systems are a few 100 days (i.e. larger than in cataclysmic 
variables) and orbital separations are lO-'^^-lO^"'^ cm. This scenario leads to more 
frequent nova explosions than in cataclysmic variables, with typical recurrence pe- 
riods smaller than 100 years. There are 10 known recurrent novae for which more 
than one outburst has been recorded and four of them are hosted by binaries with 
a red giant companion (Schaeffer 2010). 

2. THE INTERESTING CASE OF RS OPH 

The recurrent nova of the symbiotic class RS Ophiuchi had its two last erup- 
tions in 1985 and 2006. Its orbital period is 456 days, the recurrence period is 
21 years, the estimated distance is 1.6 ± 0.3 kpc, and ts = 9.5 days. 

A first issue is to find a self consistent model for the outburst of this nova, 
because it is not easy to find a combination of initial conditions leading to a 
thermonuclear runaway with such a short recurrence period. An extremely large 
initial mass and a high accretion rate are required (see Hernanz & Jose 2008). 
Initial white dwarf mass should be larger than 1.35 Mq and the mass accretion 
rate not smaller than 2 x 10^^ Mq yr~^ (see, e.g., Walder et al. 2008 and Worters 
et al. 2007, for some estimates of the expected accretion rate). 

Once the nova explodes, an expanding shock wave forms which sweeps the red 
giant wind: the system behaves as a "miniature" supernova remnant, evolving 
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Figure 2: a) Calculated postshock temperature for two values of rjinj compared 
to the RXTE and Swift data (same symbols as in Fig. 1). b) Cosmic-ray escape 
energy fraction, starting on day 6, for the best- fit acceleration efficiency Tyinj- 



much faster and being much dimmer. The characteristic timescale of evolution, 
which varies as M^/^tiwind/[£'out ('^-^/'^'•)wind], much shorter than in the super- 
nova case. Typical values for RS Oph are Moj = 3 x 10~^ Mq for the ejected 
mass, Eout — 10^^ erg for the outburst energy, and {dM/dt)v,ind = 10~^ Mq yr~^ 
and Wwind = 10 km s^^ for the red giant wind mass loss rate and terminal ve- 
locity, respectively. Typical values for core collapse supernovae of type II (SN II) 
are Mej = 10 Mq, Eout = lO^^ erg, (dM/dt)wind = lO'^ Mq yr"!, and Wwi„d = 
10 km s~^. Therefore, the characteristic time of evolution for RS Oph is about 
10~^ times that for SN II, such that the free expansion phase lasts just a few days 
in the nova remnant. 

The blast wave evolution of RS Oph during its 2006 outburst is shown in 
Figure 1 (see Tatischeff & Hernanz 2007). The time dependence of the forward 
shock velocity as deduced from IR spectroscopic observations is compared to that 
from the X-ray observations with RXTE; the usual relation for a test-particle 
strong shock Vg = [(16/3)(A:Ts)/(/imH)]^^^ was used to get shock velocities from 
X-ray fits of the postshock temperature (here k is the Boltzmann constant and 
/xmn the mean particle mass). Two questions arise from a look at Figure 1: a) 
what makes the X-ray measurements of Vg smaller than those from IR data? b) 
Why did the cooling phase started as early as 6 days, when Tg was about 10^ K 
and thus radiative cooling was not important? The answer to both question is 
cosmic rays. 

Acceleration of cosmic rays can influence the evolution of the nova remnant, 
mainly because of the energy loss due to the escape of the highest energy particles 
from the blast wave region. Then, a good agreement between the IR and X-ray 
measurements of the shock temperature can be obtained (see Figure 2a) with a 
moderate acceleration efficiency rjinj ^ lO^'* (Tatischeff & Hernanz 2007), this 
parameter being the fraction of total shocked protons in protons injected from 
the thermal pool into the diffusive acceleration process. The fraction of the total 
energy flux processed by the shock that escaped via diffusion away from the shock 
system of the highest energy particles is shown in Figure 2b. The corresponding 
energy loss rate is estimated to be ~ 2 x 10'^®(t/6 days) erg s^^ (Tatischeff & 
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Figure 3: Evolution with time from outburst of the 7-ray luminosity (in units of 
lO'^^erg s"^) corresponding to pion production and IC emission in RS Oph (2006). 

Hernanz 2007), which is approximately 100 times the bolometric luminosity of the 
postshock plasma 6 days after outburst. This means that energy loss via escape of 
accelerated particles is much more efficient than radiative losses to cool the shock. 

A prediction of the 7-ray emission associated to the accelerated particles has 
been made. The neutral pion (tt^) production has been calculated from the gas 
density in the red giant wind (i.e. {dM / dt)^ind and Uwind) and the cosmic-ray 
energy density required to explain the IR and X-ray observations. The Inverse 
Compton (IC) contribution has been estimated from the nonthermal synchrotron 
luminosity Lgyn and the ejecta luminosity Lcf Lie = ^syn x Ui-^d/ (B'^ /Sir), 
where C/rad ^cj/(47rci?^) {Rs is the shock radius and c the speed of light) 
and B is the postshock magnetic field. According to early radio detections of 
RS Oph at frequencies lower than 1.4 GHz by Kantharia et al. (2007): Lgyn ^ 
5 X 10^^{t/l day)~^-^ erg s~^. The ejecta luminosity can be assumed to be equal to 
the Eddington luminosity in novae in outburst (i.e. approximately 10^* erg s~^). 
We then find Ljc ^ ^syn and can conclude that 7-rays come mainly from tt'^ 
production (see Figure 3). 

Finally, predicted light curves for 7-ray energies > 100 MeV and > 
30 GeV are compared in Figure 4 to the Fermi/LAT sensitivities for 1 week and 1 
month observation times. The conclusion is that RS Oph (2006) would have been 
detected by Fermi/LAT. The origin of the high energy 7-rays is the collision of the 
nova blast wave with the dense wind of the red giant companion; such a situation 
can only be found in novae occurring in symbiotic binaries. 

3. V407 CYG: THE FIRST NOVA DETECTED IN HIGH ENERGY 7-RAYS 

V407 Cygni is a symbiotic Mira binary system, composed of a white dwarf and 
an asymptotic giant branch (AGB) star of the Mira type. The pulsation period 
of the Mira is approximately 750 days and the orbital period of the system is 
^43 years (Munari et al. 1990). V407 Cyg underwent an outburst on 10 March 
2010, with a spectroscopic evolution resembling that of the 2006 outburst of RS 
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Figure 4: High energy 7-ray light curves expected for RS Oph (2006), as compared 
with Fermi/LAT sensitivities for > 100 MeV. 



Oph (Shore et al. 2011). Fermi/LAT detected it on the same day and then 
repeatedly for a few days, with a peak flux in 7-rays observed between 13 to 14 
March 2010, i.e. about 3-4 days after the optical outburst (Abdo et al. 2010). 
So, our predictions for RS Oph can now be tested with real observations with 
Fermi/LAT. However, less multiwavelength observations are available for V407 
Cyg (2010) than for RS Oph (2006). 

There are two main differences between the two objects. First, V407 Cyg is 
not a standard recurrent nova, as RS Oph, since no regular eruptions were known 
before 2010 (Munari et al. 2011). Second, the orbital period is much larger for 
V407 Cyg than for RS Oph (43 years vs. 456 days); then, the separation between 
the white dwarf and the AGB star in V407 Cyg is very large (about 15 AU, 10 
times larger than in RS Oph) . This means that the shock wave needs about 7 days 
to reach a distance equal to the binary separation, and thus it should propagate 
through the red giant wind perturbed by the orbital motion. On the contrary, in 
RS Oph free expansion of the shock wave through the unperturbed wind occurred 
from day one after outburst. We also note that nonthermal radio synchrotron 
emission is not (yet) reported for V407 Cyg (2010), contrary to RS Oph (2006). 

We have made a preliminary estimate of the 7-ray flux from neutral pion decay 
in V407 Cyg (2010), getting F^{E^ > 100 MeV) ~ 10"^ photons cm'^ s-\ for an 
estimated mean postshock density for a few days after outburst of ~ 2 x 10^ cm"'^, 
a nova energy output of ~ 10'*'* erg and a distance of 2.7 kpc. This result is 
consistent with the peak 7-ray flux detected by Fermi/LAT about 3-4 days after 
the optical outburst. 

4. SUMMARY AND DISCUSSION 

Recurrent novae in symbiotic binaries are expected to accelerate particles and 
emit high-energy 7-rays detectable with Fermi/LAT, because of the shock wave 
propagation in the dense wind expelled by the red giant star. According to our 
calculations, RS Oph (2006) would have been detected by Fermi. V407 Cyg (2010) 
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was detected by Fermi, but did not behave as RS Oph regarding X-ray and radio 

emissions. Other similar systems exist in the Galaxy: eventually 1-2 symbiotic 
novae per year are expected (but not necessarily detected in the optical). 
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